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Causality constraints require a super-Hubble scale suppression in the primordial power spectrum. 
This modification is implemented and a three parameter likelihood analysis is performed of the 
COBE-DMR 4-year data with respect to the amplitude, spectral index, and suppression scale. All 
suppression length scales larger than c/Hq are consistent with the data, but scales of order 4c/ Ho are 
slightly preferred, at roughly the one-sigma level. Many non-inflation models would be consistent 
with a small suppression length scale, whereas for standard inflation models, it would require small 
e- folds. Suppression scales smaller than c/Ho are strongly excluded by the anisotropy data. 
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Causality prohibits coherence between physical phe- 
nomena with superhorizon scale separation. For density 
perturbations which are generated by local causal pro- 
cesses in a Robertson- Walker universe, the implications 
of causality have been shown to imply a suppression of 
the power spectrum which decreases faster than k for 
scales larger than the horizon Q. In non- inflationary 
cosmology the horizon and Hubble radius are about the 
same size, so that no causal mechanism could produce 
super-Hubble scale perturbations. Inflationary cosmol- 
ogy is characterized by a time period in which the hori- 
zon grew exponentially fast while the Hubble radius re- 
mained essentially constant. Thus inflation provides the 
only known causal mechanism from present theory that 
allows density perturbations to be super-Hubble scale 
with respect to the present-day Hubble radius c/Hq « 
3000 h" 1 Mpc, where h" 1 = 100/if km sec" 1 Mpc -1 . 

For large scale structure, the power spectrum of the 
primordial scalar density perturbations can be written 
as 



where V is a large rectangular volume and f{k) is the long 
wavelength suppression factor, which has been added to 
satisfy causality constraints. The modification of a sup- 
pression scale to the pure power-law power spectrum is 
mandated by causality. Under very general conditions, 
the suppression factor /(fc) can be obtained from causal- 
ity to be 
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is the wavenumber of the suppression scale. 
Causality places the the strict constraint for the suppres- 
sion index m > 4 — n. A suppression factor like eq. (|^) 
also has been found in a model with cosmic strings plus 
cold or hot dark matter Q. Eq. (|l|) contains three pa- 
rameters: the power spectrum amplitude A, the power- 
law spectral index n, and the super-Hubble suppression 
scale wavenumber k m i n . 

There is a second constraint inferable from causality, 
which is convenient from the point of view of inflation- 
ary cosmology. Models of primordial density pertur- 
bations are generally classified as either inflationary or 
non-inflationary For both types, the long wavelength 
suppression of the spectrum is a general characteristic 
which is imposed by causality, although it is typically ig- 
nored for inflation models. Since for non-inflation mod- 
els the causal horizon is about the same as the present- 
day Hubble radius, causality constraints also imply that 
k m in ~ 7r.ffo m eq. (||) . It should be noted that for 
non-inflation models, density perturbations can be both 
primordial and generated after last scattering, whereas 
eq. (|l|) strictly is valid for the former type. For many 
of the latter cases, eq. (|IJ) may still be valid based on 
||. The focus here is not on specific models, although 
this caveat is worth stating. For inflation models k m i n 
is far less constrained. For such models the largest scale 
of primordial density perturbations arose from the fluc- 
tuations which first crossed the Hubble radius during in- 
flation The factor of expansion for them is given by 
the parameter N, which is the number of e-folds of cos- 
mic scale growth. In order to solve the horizon problem, 
inflation's key purpose, the minimum required expansion 
places a lower limit of N > 50 — 70. Models of the 
standard inflationary cosmology offer no convincing rea- 
son for N to be near its lower limit and generally predict 
it to be several orders of magnitude bigger. Hence, for 
standard inflation models the super-Hubble suppression 
length scale generically is expected to be large, k m i n ~ 0. 

Previous analyses of the COBE-DMR data have tacitly 
assumed that the fundamental parameter k m i n was zero. 
In this Letter we perform a maximum likelihood fit of the 
4- year DMR data with the three parameter spectrum in 
eq. (|l|). We fix the suppression index m at various val- 
ues and determine the most likely values of the amplitude 
A, and spectral index n, in the presence of a third fun- 
damental parameter, the super-Hubble suppression scale 
wavenumber, fc m i„, and we place limits on the allowable 
range of k min . 

This analysis will consider the standard case of a flat 
universe f^o = 1 with zero cosmological constant A = 0. 
The CMBR temperature fluctuation along the direction 
unit vector n is ||] 

ST(h) _ T(h) - T _ ffg v S(k) ik . y 

T ~ T 2Vc 2 ^ fc 2 ' [ 1 
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where <5(k) is the Fourier amplitude of the density con- 
trast S(r) and y is a vector of length y = 2c/H , the 
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distance to the particle horizon in a matter dominated 
universe, which points in the direction n. The power 
spectrum is defined as 

P(fc)<W =< S(k)<Kk') >=< \S(k)\ 2 > <5 kk ,, (4) 

where the brackets denote ensemble average. P(k) will 
be taken as eqs. ([!]) and (^). 

The CMBR temperature fluctuations on the celestial 
sphere are usually expressed by spherical harmonics 

ST(h)/T = J2ai, n Y lm (h), (5) 

Int 

where Yi m (n) are the spherical harmonic functions. 
Defining a rotationally invariant coefficient Ci = 1/(2/ + 
l)E m < |a;m| 2 >, one finds from eqs. (|) 

c -=^r*^ |M%)|2 - <6) 

where ji(j£) are the spherical Bessel functions. Using eq. 
e q- © becomes 

1 2^c 4 7 1 + (k min /k) m ' 1 J 

By the above conventions, the quadrapole anisotropy is 
given as Q rms -ps = yJ^Cij where T = 2.728JT is 
the mean CMBR temperature Q . 

To simplify this initial analysis, we have not consid- 
ered temperature fluctuations produced by tensor (grav- 
itational) perturbations. Most theoretical models expect 
the CMBR anisotropics to be dominated by scalar per- 
turbations. Tensor perturbation will be subject to a long 
wavelength suppression like eq. (|2|), although the power 
spectrum index n generally should be different. 

To check for possible confusion from secondary sources 
of anisotropy, such as the Integrated Sachs- Wolfe (ISW) 
effect, we have generated fully processed power spectra 
with the code CMBFAST || for a range of the cosmo- 
logical parameters VL vac and H, for which a significant 
ISW effect is possible. The ISW effect generically in- 
creases the power in low-order multipoles and can thus 
"fill in" some of the suppression generated by a cutoff 
scale. However, this requires values of Cl vac in excess of 
~ 0.5 to be significant, and in no case can the ISW ef- 
fect mimic the effects of a small suppression length scale. 
Thus we do not further consider the ISW effect in this 
paper since it cannot decrease the significance of any pos- 
sible detection of a cutoff scale. An alternative possibility 
is that an ant i- alignment of the Galaxy quadrupole with 
the CMBR quadrupole would suppress the quadrupole 
power. However, Galaxy modeling indicates that this is 
not a significant effect |8|. 

To fit the parameters of our model power spectra, we 
will use the pixel based likelihood method (pixed based 
method) introduced in and used in COBE-DMR stud- 



ies [|10|ll • This method is predated by Gaussian likeli- 
hood fits to the 2-point angular correlation function (2- 
point method) . The primary disadvantage of the 2-point 
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method is that the 2-point correlation function is not 
Gaussian distributed, so the Gaussian likelihood fit is 
only approximate. 

In the pixel based method the covariance matrix is 
computed between map pixels i and j as 

Mij ^< >= ^ £(2Z + l)Wf QP^ ■ hj), 

i 

(8) 

where Ti is the temperature in pixel i of a map, Wj 2 is the 
experimental window function that includes the effects 
of beam smoothening and finite pixel size, Ci is given in 
eq. (Q), Pi (hi ■ hj) is the Legendre polynomial of order 
I, and hi is the unit vector towards the center of pixel 
i. For pixel temperatures that are Gaussian distributed, 
the covariance matrix fully specifies the statics of the 
temperature fluctuations. The probability of observing a 
map with pixel temperatures T, given a model C/, is 

P(T\Ci(p))dT = - " , = (9) 
V 1 lW " (2tt)^ ^ d etM(Q(p)) ^ ' 

where J is the number of pixels in the map. Assuming a 
uniform prior distribution of cosmological model param- 
eters, the probability, or likelihood function, of a given 
Ci with parameters p and a given map T is then 

-if^M-^CMp))-? 

1Hp,l "-«r' (10) 

For convenience we will denote the likelihood function 
simply as L(p). 

We have evaluated the above likelihood function using 
the model power spectra in eq. (Q) for two cases: sharp 
cutoff m = oo and minimal cutoff m = 4 — n. The results 
reported below are from the COBE correlation technique 
map, which has the best estimate of the high-latitude 
Galaxy subtracted off 0. We have tested a case with 
a map which has no residual Galaxy subtracted and the 
results are not qualitatively different. 

Full details of our tri-parameter likelihood anal- 
ysis will be reported elsewhere. Here we will 
present our results in terms of the projected like- 
lihoods L(k m in,n;Q rms -ps),L(k m m,Qrms-PS,n), and 
L(n; ps)- For a likelihood function L(p), de- 

fined by a set of parameters p = (pi,P2), the projected 
likelihood L(pi;p2) is defined as L(p) for fixed pi eval- 
uated at the most likely P2. For L(k m i n ,n;Q rrns -ps) 
the most likely (k m in , n) are given in table f with 
68% confidence level (CL) uncertainties with respect 
to the projected likelihoods L(k m i n ;n,Q rms -ps) and 
L(n;k m i n ,Q rms -ps)- Table 2 gives the projected likeli- 
hood results for both cutoff models for L{k m i n , Qrms-ps) 
under the constraint n = 1. In both tables I and 2, the 
most likely quadrapole anisotropies Qrms-ps, are given 
with 68% CL uncertainties from the likelihood evaluated 
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at the most likely values of k m i n and n (or fixed n = 1 
for table 2). These errors reflect the precision of the nor- 
malization for the specified models for a fixed shape of 
the spectrum ie. (k m im n>). Finally table 3 gives the 99% 
confidence upper limits on k m i n for unconstrained n and 
n = 1. Our results reproduce those given in |rj]jL2|1 in 
the pure power-law limit, fc m m = 0. 

Table 3 confirms a detection of a coherence length 
bigger than the Hubble radius, c/Hq, for both cutoff 
models. Tables 1 and 2 show that all values of the 
suppression length scale larger than the Hubble diam- 
eter (X m ax ~ 2c/flo) are consistent with data, although 
length scales of order 4c/ are slightly preferred. If 
the spectral index is fixed at n = 1 (Table 2), there is 
a one sigma exclusion of k m in = 0. However with n left 
unconstrained (Table 1), only the sharp cutoff model is 
found to exclude k m i n = at 68% confidence (Fig. 1). 
In no case was a two sigma exclusion of k m in = found. 
The quadrapole anisotropies Q rm .s-ps m au cases in ta- 
bles 1 and 2 are smaller than the pure power-law result 
of Qrms-ps = 15-3+2 g f*K JTT| , p~2| ] , because of the low- 
£ suppression in the spectrum. Note that our results 
for Qrms-ps are comparable to those in [Til ] where the 
quadrapole C2 was fit independent of the rest of the spec- 
trum. This suggests that the shape of the likelihood in 
our current analysis is being driven primarily by the low 
quadrupole, and the most-likely normalized spectra are 
such that the mean quadrupole in each case is compara- 
ble to the actual quadrupole in our sky. 

One can get a feeling for the relative preference be- 
tween the sharp and minimal cutoffs by examining the 
likelihood at k m i n = 0, where both model spectra are 
common. By this approach the sharp cutoff is slightly 
preferred to the minimal cutoff by a factor 1.24. Note 
that the most-likely spectral index n is closer to the scale- 
invariant limit (n = 1) with I10Il-Z6rO fernin til £1 11 it WcLS 
with k min = (n = 1.2 ± 0.3 at k min = 0,|TJ). In 
summary, a finite super-Hubble suppression length scale 
is, at best, suggested by the data. 

Two features of the power spectrum model should 
be mentioned. First, for very large suppression length 
scales, k m i n y < 1, the model spectra are virtually indis- 
tinguishable from a pure power-law spectrum. Therefore 
this analysis is most powerful at placing upper limits on 
the suppression scale wavenumber k m i n . Second, there 
is some degeneracy between k m i n and n, particularly for 
the minimal cutoff model with m — 4 — n: decreasing the 
suppression scale k m i n steepens the slope of the power 
spectrum at low spherical harmonic order £, which par- 
tially mimics the effect of increasing the spectral index 
n. Additional large and medium scale anisotropy data, 
as expected from the forthcoming satellites MAP and 
Planck, should allow us to place better limits on n so as 
to partly constrain this degeneracy. 

To cross-check the results of our likelihood analysis, 
we simulated 1000 pure power-law, scale-invariant skies, 
(k m i n ,n) = (0.0,1.0), to determined what fraction have 
likelihood functions similar to the data. For the upper 
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limit on k m i n , the Monte Carlo results firmly confirm the 
99% CL given in table 3 for both types of cutoffs. For 
the lower limit on k m i n , the Monte Carlo results indi- 
cate that a pure power-law, scale invariant universe has 
a 20% (33%) chance of spuriously imitating a universe 
with a super-Hubble suppression scale of the size favored 
by the data with a sharp (minimal) cutoff. These results 
are consistent with the likelihood analysis. In particu- 
lar they verify that the tendency found in the data for a 
small suppression length scale is not an effect of anoma- 
lously large spurious fluctuations. The Monte Carlo anal- 
ysis also indicates that increasing the signal-to-noise ratio 
with further data may give a 1-2 sigma greater discrim- 
ination of kmin = 0, but cosmic variance prohibits much 
greater significance. Thus, while the evidence for a finite 
suppression scale is not statistically significant, we find 
it suggestive enough to consider its implications. 

The standard inflation models predict k m i n ~ 0, which 
is consistent with, but not preferred by, the data. It is 
worth noting that the super-Hubble suppression behavior 
in the basic new inflation jl3| and chaotic inflation 
scenarios is closer to the sharp cutoff form. For a large 
class of non-inflationary models there is no mechanism 
for the growth of super-Hubble scale perturbations. In 
models with a "late time" cosmological phase transition 
fill , we expect k m i n y > 4 jUll. In models with topo- 
logical defects, in which the perturbations are produced 
before last-scattering, the super-Hubble suppression scale 
depends on the dynamics of the defects. For example, in 
models with cosmic strings plus hot or cold dark matter 
it was found that k m i n y ~ 2.1 — 7.9 with a cutoff be- 
havior closer to the minimal form |Q . These two types of 
non-inflation models are potentially consistent with the 
COBE data, but subhorizon evolution must be checked. 
For cosmic strings, recent simulations |l(| indicate that 
subhorizon evolution induces greater power in the low- 
order multipolcs, Cg. If further study supports this find- 
ing, it will imply that a small super-Hubble suppression 
length scale in the COBE data cannot be explained by 
cosmic strings. Both the super-Hubble suppression scale 
and cutoff should be noted in all theoretical models of 
primordial density perturbations. 

If further investigation substantiates the suggestion of 
a small super-Hubble suppression length scale, there are 
at least two possible explanations. One is that density 
perturbations produced by non-inflation models are the 
dominant contributor. The second is that density per- 
turbations are primarily produced during inflation, but 
that the number of e-folds, N, is close to its lower bound. 
In standard inflation models a small N is viewed as a fine 
tuning of the theory. Nevertheless, a small N is consis- 
tent from two other directions. Firstly, observationally 
a small N also could consistently explain a nearly-flat 
to open universe p7|| . We plan to study the suppres- 
sion scale in low density power spectrum models. Sec- 
ondly, it has been shown from standard Friedmann cos- 
mology |Q that the requirement of small N can be re- 
alized by a symmetry breaking phase transition at finite 
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temperature, during which the universe smoothly goes 
from an inflation-like stage to the radiation dominated 
stage without an intermediate reheating period. The nat- 
uralness of a small N found in this study Jul motivated 
the search in this paper for a small super-Hubble sup- 
pression length scale. Scenarios in this non-isentrephic 
inflation-like expansion regime, which was named in JlS|] 
the big-bang-like inflation regime, avoid requiring both 
localized fields on ultraflat potential surfaces and an im- 
pulsive large scale energy release during a reheating pe- 
riod Q, which must be isotropic over causally vastly- 
disconnected regimes. Both these fundamental conun- 
drums have been difficult to resolve in the standard in- 
flation picture. The theory of density perturbations ]l9| ] 
in this intermediate regime of radiation and vacuum en- 
ergy requires further study of warm inflation |20| and 
possible other mechanisms before direct comparison is 
possible with standard inflation models. The increased 
complexity of this mixed fluid regime poses several the- 
oretical questions which must be understood before firm 
predictions can be made. 

The super-Hubble suppression scale added to the 
power spectrum in this paper is required on first principle 
by any causal theory. It is more fundamental than the 
spectral index and the amplitude and much less model- 
dependent than other processes that can effect the large- 
scale anisotropy, such as the Integrated Sachs- Wolfe ef- 
fect. Moreover, the suppressing effect of causality on 
the power spectrum can produce a potentially signifi- 
cant effect on large scales that is not readily confused 
with other processes that generically boost the large-scale 
anisotropy power. This possibility calls for a reexamina- 
tion of existing parameter fits to the COBE data. 

In conclusion, we modified the primordial power spec- 
trum of density fluctuations from a pure power-law form 
to a form that includes a super-Hubble suppression scale, 
kmin, so as to properly respect causality constraints. We 
fit this spectrum to the 4- year COBE-DMR data and find 
that the data prefers a finite suppression scale, but does 
not rule out k m i n = 0. The best fit to the COBE 4-year 
data is (n,k mm y,Q rms „p S ) = (1.07±g;|§, 3.4±^, 10.9 ± 
0.7 /J.K), with a slight preference for a sharp super-Hubble 
cutoff. Upper limits on k m i n y have been firmly placed. 
We conclude that this third fundamental parameter is 
measurable from the COBE data. 

Financial support was provided in part by the U. S. 
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Sciences. 
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FIGURE CAPTIONS 

Figure 1: Relative likelihood for the suppression scale 
wavenumber, k m i n , projected over n and Qrms-PS with 
V = 2/flo- 
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Tabic 1: Maximum likelihood parameter estimates for 
unconstrained n with 68% CL uncertainties from the 
respective projected likelihood functions. 
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Table 2: Maximum likelihood parameter estimates for 
constrained n = 1 with 68% CL uncertainties from the 
respective projected likelihood functions. 
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Table 3: 99% CL upper limits on k m i n y. 
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